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We present analytical results and kinetic Monte Carlo simulations for the mobility and microscopic 
structure of solid-on-solid (SOS) interfaces driven far from equilibrium by an external force, such 
as an applied field or (electro) chemical potential difference. The interfaces evolve under a specific 
stochastic dynamic with a local energy barrier (an Arrhenius dynamic), known as the transition 
dynamics approximation (TDA). We calculate the average height of steps on the interface, the 
average interface velocity, and the skewness of the interface as functions of the driving force and 
the height of the energy barrier. We find that the microscopic interface structure depends quite 
strongly on the barrier height. As the barrier becomes higher, the local interface width decreases 
and the skewness increases, suggesting increasing short-range correlations between the step heights. 

PACS numbers: 68.35.Ct 75.60. Jk 68.43. Hn 05.10.Ln 



I. INTRODUCTION 

A large fraction of the processes essential to some of the most important technologies in our society involve the 
physical and chemical properties of interfaces: an interface is the place where molecules from different phases come 
into contact, and where chemical bonds are formed. Catalytic reactions and semiconductor devices are only two 
examples of systems where the role of surfaces and interfaces is fundamental P, 0j- It is therefore of the utmost 
importance to understand the fundamental processes occurring at surfaces and interfaces in order to design more 
efficient mechanisms to contribute to the technological progress. 

Since the specific dynamics of the evolution of surfaces and interfaces are in general not known, it is useful to 
construct stochastic models that reproduce the essential features of the physical system. Until recently it was com- 
monly believed that dynamics that respect the same conservation laws and all obey detailed balance, give essentially 
the same qualitative behavior. However, recent results clearly show that there are nonconservative dynamics that, 
although they all obey detailed balance and the same conservation laws, lead to very different interface microstruc- 
tures [H S H S ■ Since many interface properties, such as mobility and chemical activity, are determined by the 
microstructure, great care must therefore be taken in selecting the appropriate dynamic for the physical or chemical 
system of interest. 

In this work we present analytical results and kinetic Monte Carlo (MC) simulations for the mobility and microscopic 
structure of (l-l-I)-dimensional solid-on-solid (SOS) interfaces 8] in an Ising lattice-gas model, which are moving far 
from equilibrium under an external force, such as an applied magnetic or electric field or an (electro) chemical potential 
difference. The system evolves under an Arrhenius dynamic which includes a microscopic barrier that represents a 
transition state inserted between the states allowed in the Hamiltonian. Arrhenius dynamics are appropriate when 
discrete Ising or lattice-gas models are used to simulate dynamics in an underlying, continuous potential. In such cases 
the transition state in the Arrhenius dynamic could represent a saddle point in a corrugation potential for particle 
diffusion, or a high-energy state associated with a transitional s pin state [3.lfC|. The dynamic considered in this paper 
is the two-step transition dynamic approximation (TDA) |lOl If f|. The transition-state energy is approximated by 
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where Ef and -B, are the initial and final energies, and A is the microscopic energy barrier. In electrochemical 
applications, this corresponds to the symmetric Butler- Volmer approximation |l2j |. The transition rates for the TDA 
are given by 

Wtda = (2) 

1 + cxp[/3(£ T - Ei)\ 1 + exp[(3(E f - E T )\ 1 ; 

The TDA belongs to the class of dynamics defined as "hard," in which the transition rates cannot be factorized into 
one part that depends only on the interaction energy and another that depends only on the field energy. Dynamics 
that have this factorization property are defined as "soft" When app lied to kinetic Ising lattice-gas models, 

soft and hard Arrhenius dynamics give very different nucleation rates . It has also recently been shown that 
SOS interfaces evolving under soft Arrhenius dynamics have quite different structures that those evolving under hard 
Arrhenius dynamics 7J. In soft Arrhenius dynamics, the barrier height only results in a change in the overall time 
scale of the simulation and has no effect on the interface microstructure Q • In the present work we focus on the effect 
of the barrier height A on the mobility and microstructure of interfaces moving under the hard TDA dynamic. 



II. THE MODEL 



In this section we give a brief description of the SOS model and the analytical theory for the interface microstructure 
and mobility. Details can be found in Refs. 0,0, 0,0- 

The SOS interfaces are described by the nearest-neighbor Ising Hamiltonian, 

7~l = — S x<y (J x S x +l,y + JyS x ,y+l + H) , (3) 

where s x ^ y — ±1, ^ x runs over all sites, and the applied field H is the driving force. The interface is introduced by 
fixing s XtV = +1 and —1 for large negative and positive y, respectively. We take H > 0, such that the interface on 
average moves in the positive y direction. This Ising model is equivalent to a lattice-gas model with local occupation 
variables c XtV E {0, 1} (for the exact relations between the parameters in the two models, see, e.g., Ref. 7J). 

The interface is described by a single-valued integer function h(x) with steps 8{x) = h(x + 1/2) — h[x — 1/2) at 
integer values of x. An example of an SOS interface is shown in Fig. 1. The spins in the anisotropic square-lattice 
Ising model can be divided into classes, labeled jks, where j and k are the number of broken bonds between the spin 
and its nearest neighbors in the x and y directions, respectively, and s is the spin value. Only spins that belong to 
the classes 01±, 11±, and 21± can flip in the SOS interface. Whenever a spin flips from —1 to +1, the corresponding 
column of the interface advances by one lattice constant in the y direction. Conversely, the column recedes by one 
lattice constant when a spin flips from —1 to +1. In this approximation the spin-class populations on both sides of 
the interface are equal, and the contribution to the mean velocity in the y direction from sites in the classes jk— and 
jk+ becomes 

(v y (jk)) = W (/3SE(jk-),/3A) - W (/38E(jk+),/3A) . (4) 

where SE(jk±) is the change in energy due to flipping a spin in the class jk±, and the transition rates W are 
calculated from Eq. (2). The mean propagation velocity perpendicular to the interface is defined as 

(v ± (T,H,<j>)) =coscf>J2n3ks))(v y (jk)) , (5) 

where <f> is the overall angle between the interface and the x axis. In this mean-field approximation, the heights of 
the individual steps are assumed to be statistically independent. As a consequence, the mean spin-class populations 
(n(jfcs)) are obtained from the product of the step-height probability density function (pdf) for 8(x) and S(x + 1). 
The pdf is given by the interaction energy corresponding to the \S(x)\ broken J x bonds between spins in the columns 
centered at (x — 1/2) and (x + 1/2) as 

p[5(x)} = Zi^X^^ e 7W)<5(a:) . (6) 

The factor X determines the width of the pdf; in equilibrium it is simply the Boltzmann factor, e~ 2/3Jn: . The quantity 
7(0) is a Lagrange multiplier that maintains the mean step height at an x-independent value, (S(x)) — tan</>. The 
partition function is 

o= — 00 
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where 7(0) is given by 
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A non-linear response approximation based on a dynamic mean-field approximation for the equation of motion of the 
single-step pdf, together with a detailed-balance argument for the stationary state, gives the following expression for 
the field-dependent X{T,H) dH, 

X(T H) e-^ h- 2fm W(2l-) + e^W(21+) \ 1/2 

X ^ H )- e { W {2l-) + W{2l+) / ' (9) 

where the W(jks) are the single-site transition rates associated with the flipping of a spin of the class jks. For the 
details of the calculation see Refs. HQ. For the TDA dynamic defined in Eq. (2) one gets 

y (T m _ -20 Jx \ ^ cosh(2/3g) + e~ 2 ^ + 2 cosh(/3A) costing) 1 1/2 
TDAl ' ' ' \ e -W» cosh(2/?#) + e 2 '"* + 2 cosh(/3A) cosh(/3H) J ' 1 j 



III. RESULTS 

The numerical simulations were performed with a continuous-time n-fold way rejection-free algorithm 14] at T = 
0.2T C (where T c w 2. 269 J is the critical temperature for the isotropic, square-lattice Ising model) with L x = 10 000 
and isotropic interactions, J x — J y = J. In Fig. 2 we show the results for the average stationary step height (\5\) vs 
H, for different barrier values. The agreement between the analytical results and the MC data is very good for the 
smaller values of A. We believe that the differences between the analytical and the MC results for the higher values 
of A, particularly for small H, may be related to the fact that the assumption of independent step heights is not valid 
for larger barriers due to the increasing presence of short-range correlations. Figure 2 indicates that the local width 
of the interface depends strongly on H and on the height of the barrier. The larger the barrier, the smaller the local 
interface width. We also note that the agreement is not very good for H < A. However, the agreement is almost 
exact for H — A, and it remains quite reasonable for H > A. 

The dependence of the perpendicular interface velocity on the field is shown in Fig. 3 for several values of A. The 
agreement between the analytical and MC results for the velocity appears to be very good for all values of A and H . 
This is, however, somewhat of an optical illusion. The good agreement occurs for the higher velocities, which all are 
attained for H > A. If one concentrates on the low velocities attained for H < A, one indeed finds that the theory 
significantly underestimates the simulation data. This is as expected, since the dominant contribution to the velocity 
at this low temperature is from spins in class lis, whose abundance in the simulated interface is larger than expected 
from the theory. As predicted by the theory, the interface velocity of systems that evolve under the TDA is bounded 
by unity, contrary to the case of the soft Arrhenius dynamic considered in Ref . Q , for which it increases exponentially 
with H. From the figure it is clear that, as the barrier increases, the velocity increases more slowly with H, and it 
reaches its maximum value at higher fields. 

As mentioned above, the analytical predictions for the class populations are based on the mean-field assumption 
that the steps are statistically independent. Within this approximation, the average populations of each class of spins 
must be the same in front of and behind the interface. However, the simulation results clearly show that this is not 
the case in general. As an example, the six mean class populations -(n(01s)}, (n(lls)}, and (n(21s)) with s = ±1 - 
are shown vs H in Fig. 4, for A = 2 J. Again we notice that, for this relative high value of A, there is not a complete 
agreement between the theoretical and the MC data at intermediate fields. The results indicate the presence of 
short-range correlations between neighboring steps that cause differences between the spin populations on the leading 
and trailing edges of the interface. When (n.(21— )) > (n(21+)), the interface is characterized by a broadening of 
protrusions on the leading edge ( "hilltops" ) . The relative skewness can be quantified by the following quantity 

_ (n(21-))-(n(21+)) 
P (n(21-)) + (n(21+)) ' 1 ) 

This skewness parameter is shown in Fig. 5 for different values of the barrier height. Notice that the relative skewness 
shows a dramatic increase as the barrier height becomes larger, particularly for low and intermediate fields. These 
results suggest that the assumption of independent step heights fails for systems with high barriers. Skewness has 
also been observed in other SOS-type models [lallrill7|. 
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IV. CONCLUSIONS 



In this work we have presented numerical and analytical results for the mobility and interface microstructure of an 
SOS interface, which is driven far from equilibrium by an applied field. The system evolves according to an Arrhcnius 
dynamic known as the transition dynamics approximation (TDA). The Arrhcnius dynamic interposes between the 
Ising lattice-gas states a transition state representing a local energy barrier. In this study we focus on the effects of 
the barrier height on the interface microstructure and mobility. The analytical results are obtained with a non-linear 
response theory, which assumes that there are no correlations between the heights of different steps on the interface. 

We found that the microscopic properties of the interface depend quite strongly on the barrier height. The interfaces 
are less wide, move more slowly, and are more skewed as the barrier height increases. The agreement between 
theory and simulation is very good for the higher interface velocities (Fig. 3). For low barriers, the local interface 
width is also in good agreement, but for larger A, the agreement only becomes good for H > A (Fig. 2). For 
higher barriers, the results suggest that short-range correlations become important (Fig. 4, 5). To account for such 
correlation effects, an improved theory would be needed. Finally we emphasize that the disagreements between theory 
and simulations presented in this paper represent a "worst-case scenario." The agreement improves considerably for 
higher temperatures (below T c ), and at T = 0.6Tc it is quite good for all barrier heights and fields studied. 
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FIG. 1: A segment of an SOS interface y = h(x) between a positively magnetized phase (or "solid" phase in the lattice-gas 
picture) below and a negative (or "fluid") phase above. Sites in the uniform bulk phases are 00— and 00+. From Ref. U 
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FIG. 2: Average stationary step height (\5\) vs H for 0=0 at T = 0.2T C . The curves represent the analytical results. The MC 
data were obtained directly by summation over the simulated single-step pdfs. 
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FIG. 3: The average stationary normal interface velocity (v±) vs H at T = 0.2T C , calculated for several values of the barrier 
A. The symbols represent MC data, and the solid curves analytical results. 




FIG. 4: Mean stationary class populations {n(jks)} vs H/J at T = 0.2T C , for A = 2 J. The symbols represent MC data, and 
the solid curves analytical results. 
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FIG. 5: MC data for the relative skewness p, defined in Eq. (6), vs H at T — 0.2T C , for different values of the barrier height A. 
The solid lines are merely guides to the eye. 
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